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OA , We calculate thermodynamic quantities and transport coefficients of a hadronic gas using 

<—, ■ a Monte Carlo hadronic collision event generator URASiMA. The obtained shear viscosity 

to entropy density ratio of the hadronic gas is approximately 0.3 - 0.6 in the range T = 100 
h— j ' - 180MeV and is quite insensitive to the baryon number density. 

\Q ■ 

§1. Introduction 

Hydrodynamical models are among the most thoroughly established models of 
multiple production phenomenal In particular for RHIC experiments, hydrody- 
q . namical models are accepted as the most successful models to reproduce V2 data.® 

To this time, in almost all hydrodynamical models, viscosities and heat conductivity 
have been ignored for simplicity. Teaney® attempted to estimate the size of the effect 
that may appear in the present hydrodynamical models applied to RHIC experiment 
in the case that the shear viscosity is taken into account, and he concluded that the 
viscous term in the Navier-Stokes equation must be smaller than unity in order to 
retain good agreement between the blast wave model with viscous corrections and 
experimental results. A peculiarly small ratio of the viscosity to entropy has been 
identified as a possible special feature of the QGP liquid.™ Here we evaluate the 
J> . shear viscosity and entropy for a hadronic gas. 

| Calculating macroscopic material constants on the basis of a microscopic picture 

is one of the most important tasks in theoretical physics. This is, however, very 



difficult in general. Particularly in the case of strongly interacting hadrons, there 
exists no perturbation theory and no systematic method to calculate the relevant 
quantities. 

In this letter, we evaluate the transport coefficients of dense, hot hadronic matter 
employing a hadro-molecular dynamic calculation that an event generator URASiMA 
(Ultra-Relativistic AA collision Simulator based on Multiple Scattering Algorithm) .® 
The obtained heat conductivity and shear viscosity exhibit a temperature depen- 
dence of ~ T 5 , which is stronger than the result obtained from naive dimensional 
analysis, ~ T 3 . Muronga attempted a similar calculation for a pure meson gas using 
UrQMD, in which case, several interactions must be switched off by hand in order 
to realize the stationary state.® In the present letter, we investigate a finite baryon 
number system with all hadronic interactions, including those of baryons. 
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Table I. Baryons, mesons and their resonances included in the URASiMA. 



nucleon 


N93S 


Al440 


N152O 




Al650 


^1675 


A168O iVj.720 


A 




^1600 


Zil620 


^1700 


^1905 


^1910 


^1950 


meson 


TV 


T]s 


C"800 


P770 









§2. Hadro-molecular dynamics 

In order to construct a statistical ensemble state, we carry out a molecular 
dynamical calculation of hadrons in a box. The time evolution of the system is 
simulated with the Monte-Carlo collision event generator URASiMA. One of the 
present authors (N. S.) improved and tuned the parameters of URASiMA so as to 
realize an equilibrium state of hadrons for temperatures in the rande 100 - 180 MeV 
for nuclear densities near twice the normal nuclear density P In Ref. 7), equilibration 
of the system is discussed in detail and thermodynamical quantities are analyzed. 
Several kinds of diffusion constants have already investigated on the basis of the 
ensembles 

The present version of the URASiMA contains the particles listed in the Table 
I. This version is the so-called "two-flavor" and "low-energy" version, which ignores 
anti-baryons and strangeness. As a hadronic gas model, this model covers the region 
characterized by temperature lower than Kaon's mass and chemical potential lower 
than the baryon mass. The region on which we focus our consideration is below T c 
(~ 200 MeV) and up to twice the normal nuclear baryon number density. We believe 
that this region is well within the regime in which the present version of URASiMA 
is valid. 

Thermodynamical quantities of the systems of interest are investigated in Ref. 7) . 
Figure 1 displays the obtained equation of state (the pressure P as a function of the 
energy density e). Throughout our simulation, the box size is fixed as 10 3 fm 3 . 




Fig. 1. Equation of state of a hadronic gas Fi S- 2 - Entropy of a hadronic gas obtained 
obtained using URASiMA. using URASiMA. 
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The black circles represent the normal nuclear density ni,o (n^o = 0.157 (1/fm 3 )), 
the white boxes represent 1.5 times n^o, and the triangles represent 2 times n^o- We 
note from Fig. 1 that the bulk viscosity of the hadronic gas almost vanishes. Because 
the current P' which corresponds to the bulk viscosity is given by^ 

OP 

P'=p- B ^ (1) 

P' vanishes if the pressure is proportional to energy density. 

Figure 2 displays the entropy density calculated in Ref. 7) as S = ^ phi/9, where 
p is the ensemble averaged density function and ^ represents the sum over phase 
space, respectively. Throughout our simulation, the baryon number is conserved. 
The entropy density grows more rapidly than that in the free hadron model (the 
curves in the figure) for the temperatures approximately equal to or larger then the 
pion mass. 

§3. Relaxation of the currents 

According to linear response theory, a transport coefficient is obtained from the 
correlation of the corresponding current.^ The shear viscosity and heat conductivity 
are obtained from the correlation of the heat current and stress tensor Hij as^ 

n= hm 1 / d 3 x' [ dt'e-^-^iP^x)^^',^)), (2) 
r), = „lim o | [d 3 x'[ dt'e-^-^^y^x),^^,^)), (3) 



e- 



-oo 



where Pi and 7Tjj are the spatial-temporal and spatial-spatial components of the 
energy-momentum tensor, respectively. 

In the hadro-molecular dynamics described by URASiMA, the energy-momentum 
tensor is given by the sum of that of each particle: 



T^{x,t) = Y J TK5{x-x {l) {t)) (4) 



Here, the subscript I indexes the particles. Average is taken as an ensemble average: 

^ ^ number of ensembles ^ ^ ^ 

ensemble 



Correlation is assumed to exist only between the same particle in a state. 



*) 



*' A detailed explanation of the procedure will be given in Ref. 16). 
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n„= 0.157 [tan 



. 10 fm 

B = 0.231 [frrf 3 ] 



o" 

■ 

• o 



Figures 3 displays the correlation 
function of the stress tensor 7Tjj , which 
represents the relaxation of the vis- 
cous shear tensor and appears in the 
integrand of Eq. (j3J). At temperatures 
near the pion mass, mesonic degrees 
of freedom gradually come to domi- 
nate transport phenomena. Detailed 
analyses of the separate contributions 
of the baryonic sector and mesonic 
sector will be reported elsewhere.^ 

In each case, for values of t 
greater than about t = 20 fm, the 
correlation function becomes smaller 
than ~ 10 -3 and can no longer be 
distinguished from the noise. Because 
we are interested in phenomena on the 
hadronic scale, in actual calculations, 
we cut off the infinite integration in 
Eqs. (2) and (3) at t = 20 fm by hand. 

It is known that the final results are not sensitive to this cut off value. 

§4. Transport coefficients 



IE = 0.370 GeV/fm 
) e m 0.463 GeV/fm s 
1 6 = 0.556 GeV/fm s 
] 6 = 0.648 GeV/fm 3 
>£. 0.741 GeV/fm s 
>6 = 0.833 GeV/fm s 

n„= 0.315 [taf a ] 



' O 

i ♦ o 
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Fig. 3. Correlation function of the viscous shear 
tensor tt x ,v 



Figure 4 displays values of the calculated shear viscosity. Though the temper- 
ature range is not large, it is seen that the temperature dependence of r\ is ~ T 5 , 
which is significantly stronger than that obtained from naive dimensional analy- 




Fig. 4. Shear viscosity r\ a as a function of temperature. 
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• n B = 0.157 fm 
3 3 ■ n B = 0.231 fm" 3 

V = 10 fm An B = 0.315fm" 3 

-i — ■ — i — ■ — i 1.5 i . , . , . , ^- 




T [MeV] In T 

Fig. 5. Heat conductivity function of temperature. 



sis, ~ T 3 . This strong temperature dependence is caused by the liberation of the 
pion degrees of freedom, because the temperature is slightly above the pion mass. 
A similar temperature dependence also appears in the heat conductivity shown in 
Fig. 5. The somewhat peculiar behavior of the heat conductivity in the region of very 
low temperature can be attributed to the contribution of the baryons, whose mass 
is larger than the temperature and whose number must be conserved. A detailed 
investigation will be presented elsewhere.^ 

The obtained shear viscosity is approximately 0.2 - 0.6 [GeV/fm 2 ] and is al- 
most independent of the baryon number density. Though the system treated by the 
URASiMA is not a simple pion gas, the transport of stress is dominated by pions, 
because the temperature of the system is sufficiently smaller than the masses of 
hadrons other than pions. 

About twenty years ago, Gavin calculated the shear viscosity of a pion gaa^ 
and obtained 0.05 - 0.15 [GeV/fm 2 ], which is approximately one-fourth as large as 
our result. According to the recent calculation by Dobado and Llanes-Estradaf^ 
the shear viscosity of a pion gas is approximately 0.4 [GeV 3 ]( ~ 10 [GeV/fm 2 ]) if a 
constant scattering amplitude is assumed (see Fig. 2 of Ref. 13)) and 0.004 [GeV 3 ]( ~ 
0.1 [GeV/fm 2 ]) in the simple analytical phase shift model (see Fig. 3 of Ref. 13)). Our 
model gives a result between these values. For a pure meson gas, using UrQMD in 
which several interactions are modified to equilibrate the system, Muronga obtained 
a shear viscosity which is almost one-half of our result 

§5. r) s to entropy density ratio 

Usually, the Reynolds number R is evaluated to estimate how large an effect is 
caused by viscous terms. If Bjorken's 1+1 dimensional Scaling solution is adopted 
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as a particular solution, the inverse Reynolds number, R 1 , is given by^ 



R 



-i 



Vs + 3^ 

r(e + P) 



(7) 



s rT' 



(8) 




where r is the proper time. 

Instead of the inverse Reynolds num- 
ber, which depends on the solution of hy- 
drodynamics, we can use the shear viscos- 
ity to entropy density ratio for simplisityP^ 
The shear viscosity to entropy density ra- 
tio, rjs/s, obtained in our simulation is plot- 
ted in Fig. 6. It is seen that rj s /s is ap- 
proximately 0.4 and quite insensitive to the 
baryon number density. The peculiarly 
small value rj s /s for QCD matter has been 
pointed out by several authors, 41 but even 
in the case of hadronic gas, the value of T] s j s 
obtained in our simulation is approximately 
order unity or smaller. The URASiMA con- 
sists of ordinary hadronic collisions onlyP 
that is, there exists no special strong corre- 
lation, such as strings and potentials. The 

main difference between URASiMA and ordinary non-relativistic molecular dynamic 
simulations lies in the relativistic properties of the system. Therefore, we believe that 
the result rj s /s < 1 may be a common feature of a relativistic gas. More detailed 
analysis is now in progress P^* 
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Fig. 6. Viscosity to entropy density ratio 
of a hadronic gas. The black circles 
represent the normal nuclear density 
Tito, the boxes represent 1.5 times nio, 
and the triangles represent 2 times n^o • 
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